.t  NT  Of  j~ 


I  run ' 


# 


Ov  o 

^/v  s  pO 


AFTERBODY  DRAG 

VOLUME  2  -  JET  INTERFERENCE  EFFECTS 
ON  SUBSONIC  BOATTAIL  DRAG  A 


Peter  R.  Payne 
Payne,  Incorporated 


rs 

K,s> 


APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED 


AVIATION  AND  SURFACE  EFFECTS  DEPARTMENT 


DTNSRDC/ASED-80/12 


June  1980 


IIM'JII 


TAYLOR 
NAVAL 
SHIP 
RESEARCH 


umai  lUu  ai 


BUllli 


1 1 1 1 


La 

Kr2 


BE  TH  ESDA 
MARYLAND 
20084 


80 


I  0-10 

^  ri  _  . 


AECUfllTV  CLASSIFICATION  OF  THIS  PAGE  flFhwl  Dele  Entered) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


I.  RECIPIENT’S  CATALOO  NUMBER 


a  ijf/lTi 


r0T: 


t  cm. 


TYPE  OF  REPORT  *  PERIOD  COVEREO 


VOLUME  2 i  JET  INTERFERENCE  EFFECTS  ON 
"Subsonic  BOTTAfL  drag  ,  < 


m 


t.  PERFORMING  ORGANIZATION  NAME  AND  AOOREIS 


^ y  Final t 


mmm 


)  N00014- 7^0039  j 


AR?A  VwoVk'unFtVuMBER^-  ™K 

(See  reverse  side) 


Payne  Incorporated 
Annapolis,  Maryland  21401 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  ,  j  A 

David  W.  Taylor  Naval  Ship  R&D  Center  •  // 

Aviation  and  Surface  Effects  Department 
Bethesda,  Maryland  20084 


MONITORING  AGENCY  NAME  «  ADDRESS/-//  dltloronUn^.  QonltlWoM.OHIco)  IS.  SECURITY  CLASS,  (of  Oil*  report; 

Naval  Air  Systems  Command  / >  o  )  ] 

Washington  D.C.  20361  \!'C/  /A  !  UNCLASSIFIED 

Naval  Weapons  Center  >  f  ^  i  is#,  declassification/ down  grading 

China  Lake,  California  93555  ' _  f - - - - - 


IS.  DISTRIBUTION  STATEMENT  (ot  thto  Ropotl)  '  /£  j  7  ^  If-  A S>  S>  &(  A 


IS.  NUMBER  OF  PAGES 

75 


APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED 


St  <?t 


IS.  KEY  WOROS  (Continue  on  rororoo  e/de  it  neeeelerr  end  Identity  by  block  number; 

Base  Pressure 
Boattail  Drag 
Jet  Interference 
Subsonic  Flow 


20.  ABSTRACT  (Continue  on  reverie  e/de  II  neceeearp  end  Identity  by  block  number; 

3  The  results  of  the  afterbody  drag  study  are  presented  in  four  volumes, — 
Volume  1:  Drag  of  Conical  and  Circular  Arc  Afterbodies;  Volume  2:  Jet 
Interference  Effects  on  Subsonic  Boattail  Drag;  Volume  3:  Literature  Survey; 
and  Volume  4:  Data  and  Analysis 

^  Volume  2  includes  a  method  of  analysis  as  a  base  from  which  a  more 
detailed  analysis  of  afterbody  drag  can  be  developed.  _ 


DD  , 1473 


EDITION  OF  I  NOV  ••  IS  OBSOLETE 
S/N  0102-LF-01  4-6601 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wien  1 


3  yy  yy/ 


UNCLASSIFIED _ 

SKCUNITY  CLASSIFICATION  OF  'THIS  PAOB  <**«■  D«< 


(Block  10) 

Program  Elements: 

63361N 

62332N 

62241M 

Task  Area  Numbers: 
W15X20000 
F32.322.203 
WF41.421.201 

Work  Unit  Numbers: 
1660-234 
1660-235 


_ UNCLASSIFIED _ 

SKCUNITV  CLASSIFICATION  OFTMIS  FAOKOFNwi  DM*  Xnl#rA« 


TABLE  OF  CONTENTS 


Page 


ABSTRACT .  1 

ADMINISTRATIVE  INFORMATION .  1 

INTRODUCTION .  2 

ENTRAINMENT  DRAG  ON  THE  BOATTAIL .  5 

PLUME  INTERFERENCE .  25 

SUMMARY  AND  CONCLUSIONS .  44 

REFERENCES .  45 

APPENDIX  A  -  DERIVATION  OF  THE  SUBSONIC  SLENDER  BODY 

EQUATIONS  FOR  A  BODY  AT  ZERO  INCIDENCE .  47 

APPENDIX  B  -  PROGRAM  "PLUME3" .  63 


Accession  For 

ntis  Gam 

DDC  TAB 

Unannounced 

Justification 


p.v 

Dist.ribut1.cn/  _  ... 

Avni1 

ahU'tv  C-'  ,r  S 

Dist 

fi 

Avail  end /or 
special 

LIST  OF  SYMBOLS 


tan  0 

Sonic  throat  area 

I  2 

Static  pressure  coefficient  =  Ap/-^^ 

Velocity  mixing  coefficient  defined  in  equation  (4) 
Entrainment  coefficient  defined  in  equation  (5) 

Free  stream  Mach  number  uQ/a0 

Mach  number  of  the  plume  boundary ,  equation  (24) 

Entrainment  ratio  defined  in  equation  (5) 

Local  static  pressure 

Static  pressure  at  the  jet  exit,  equations  (26)  and  (27) 
Total  pressure  in  the  jet 
Free-stream  static  pressure 
P  -  Poo 

Strength  of  an  axial  source  or  sink  line  per  unit  length 
A  jet  plume  curvature  parameter,  defined  in  Figure  19 
Radial  distance  in  cylindrical  coordinates 
Potential  "core"  radius  ,  defined  after  equation  (4) 

Jet  radius  (see  Figure  3) 

Maximum  body  radius  (see  Figure  3) 

Jet  velocity.  Generally,  u^  =  u^Q 
Jet  "core"  velocity,  assumed  uniform 
Free-stream  velocity 

Change  in  fluid  velocity  parallel  to  the  flow  axis 

Radial  velocity  in  cylindrical  coordinates 

Change  in  radial  fluid  velocity 

Distance  along  the  axis  of  symmetry 

Downstream  distance  from  the  jet  exit  (see  Figure  3) 

Location  of  maximum  plume  diameter 

A  dummy  variable 

A  transformation  variable  with  various  meanings  which  are 
defined  where  it  is  employed 

Nozzle  angle  ,  defined  in  Figure  19 


iv 


LIST  OF  SYMBOLS  (continued) 


Y  Ratio  of  specific  heats 

6j  Initial  turning  angle  of  jet,  defined  in  Figure  19 

e  Virtual  kinematic  viscosity 

fl  Surface  angle  at  the  jet  end  of  the  boattail  (see  Figure  3) 

*  (uj/u0)  ‘  1 

v  (k/2)X 

v  Prandtl-Meyer  expansion  angle  defined  in  equations  (22)  and  (23) 

P  Fluid  mass  density 

♦  The  "shape  factor"  defined  in  equations  (17)  and  (18) 

i|>  gA  =  g  tan  0 


Suffixes 

BE  On  the  boattail,  due  to  entrainment 

gP  On  the  boattail,  due  to  jet  pluming 

o  Free  stream 

j  Jet 

t  Total 


v 


ABSTRACT 


The  results  of  the  afterbody  drag  study  are  presented  in 
four  volumes  —  Volume  1:  Drag  of  Conical  and  Circular  Arc 
Afterbodies;  Volume  2:  Jet  Interference  Effects  on  Subsonic 
Boat tail  Drag;  Volume  3:  Literature  Survey;  and  Volume  4: 

Data  and  Analysis 

Volume  2  includes  a  method  of  analysis  as  a  base  from  which 
a  more  detailed  analysis  of  afterbody  drag  can  be  developed. 

ADMINISTRATIVE  INFORMATION 

The  work  reported  was  performed  for  the  David  W.  Taylor  Naval  Ship 
Research  and  Development  Center  by  Payne ,  Incorporated  under  ONR 
Contract  N0014-77-C-0039  as  part  of  an  evaluation  of  afterbody  drag.  The 
afterbody  drag  project  was  supported  by  the  Naval  Air  Systems  Command  and 
the  Naval  Weapons  Center  under  Program  Elements  63361N,  62332N,  and 
62241N ;  Task  Areas  W15X20000,  F32.322.203,  and  WF41.421.201;  and  Work 
Units  1660-234  and  1660-235. 

The  views  and  conclusions  contained  in  this  document  are  not  necessarily 
the  official  policies  either  expressed  or  implied,  of  the  David  W.  Taylor 
Naval  Ship  Research  and  Development  Center  or  the  U.S.  Government.  The 
format  of  this  report  is  that  of  Payne,  Incorporated. 


INTRODUCTION 


This  report  is  concerned  with  the  effect  of  a  propulsive  jet  upon 
the  static  pressure  distribution  over  the  boattail  from  which  it  issues, 
and  hence  its  effect  on  the  boattail  drag.  (We  do  not  consider  "base  drag" 
and  formally  exclude  it  by  assuming  zero  base  area.)  It  is  a  fundamental 
analysis  in  the  sense  that  there  are  no  "factors  determined  from  experiment" 
in  the  usual  sense,  but  on  the  other  hand,  no  claims  to  comprehensiveness 

are  made.  The  work  presented  is  offered  as  a  basis  from  which  a  more  com¬ 
prehensive  analysis  can  be  developed,  in  that  it  does  quantify  a  reasonably 

meaningful  physical  picture.  The  theory  does  not  apply  if  the  flow  is 

separated. 

Figure  1  depicts  a  typical  variation  of  CD  with  the  jet  to  free 
stream  velocity  ratio,  which  is  a  function  of  ® 

Nozzle  pressure  ratio  (NPR) 

Jet  temperature 

Ratio  of  specific  heats  in  the  jet  (y) . 

Figure  2  gives  the  relationship  between  jet  Mach  number  Mj  and  NPR. 

Referring  to  Figure  1,  we  see  that  the  jet-off  drag  is  typically 
high  and  falls  rapidly  when  quite  small  quantities  of  fluid  are  allowed  to 
flow  from  the  nozzle.  This  region  has  been  explored  quite  extensively, 
from  an  experimental  point  of  view,  although  little  theoretical  analysis  has 
been  undertaken.  In  this  region,  the  flow  field  is  so  different  from  that 
pertaining  when  the  propulsive  jet  is  present  that  there  seems  little  hope 
of  correlating  the  two  regimes,  as  some  writers  have  attempted. 

The  region  between  "base  bleed"  and  u.  =  u  has  been  little  ex- 

j  0 

plored,  either  theoretically  or  experimentally,  probably  because  it  is  of 
little  practical  importance.  When  the  jet  velocity  equals  the  free  stream 
velocity,  we  have  the  "true"  boattail  drag  without  interference.  (This  is 
approximation  because  the  static  pressure  behind  the  nozzle  is  locally 
a  little  above  ambient,  so  the  velocity  is  therefore  a  little  below  free 
stream.)  Again,  there  is  very  little  experimental  data  available  for  this 
region;  perhaps  because  the  experimentalists  haven't  fully  appreciated  its 
significance  as  a  "benchmark"  data  point. 
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NOZZLE  PRESSURE  RATIO  (NPR) 

igure  2.  The  Variation  Between  Jet  Mach  Number  (M.)  and  Nozzle  Pressure  Ratio  (NPR)  When 
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At  higher  jet  velocities,  free  air  is  entrained  into  the  jet, 
which  therefore  acts  as  a  "sink,"  and  increases  the  velocity  over  the  boat- 
tail,  and  therefore  its  drag.  The  theory  of  this  report  suggests  that  this 
is  a  linear  effect  and  that  the  boattail  drag  will  increase  linearly  with 
jet  velocity  until  it  reaches  sonic  velocity.  This  does  seem  to  be  borne 
out  by  experiment;  but  in  most  published  studies,  it's  rare  to  find  even 
three  experimental  points  in  this  region  for  a  given  configuration. 

As  the  jet  passes  through  sonic  speed,  a  number  of  second  order 
phenomena  can  cause  discrepancies  between  one  test  configuration  and  another. 
The  next  major  trend,  however,  is  caused  by  the  exhaust  pluming,  as  its  speed 
is  further  increased.  This  has  the  same  effect  as  a  distributed  source  in 
the  wake,  and  reduces  the  flow  velocity  over  the  boattail,  leading  to  an 
increase  in  pressure  and  a  reduction  in  drag.  So  the  plume  effect  starts 
to  cancel  out  the  drag  increment  due  to  entrainment.  Eventually,  the  plume 
is  large  enough  for  the  boattail  drag  to  become  negative. 

The  effects  are  most  marked  when  there  is  no  base  area,  and  this 
is  the  case  we  have  considered.  A  base  tends  to  insulate  the  boattail 
from  these  effects,  and  they  diminish  with  increasing  base  area.  Further 
study  is  needed  for  the  more  complex  case  of  a  finite  base  with  a  jet 
issuing  from  its  center. 

ENTRAINMENT  DRAG  ON  THE  BOATTAIL 

BOUNDARY  LAYER  DISPLACEMENT  SURFACE 


Figure  3.  The  Effect  of  Propulsive  Jet  Entrainment  on  the 
Boundary  Layer  Displacement  Surface  (M^  -  0.9). 


Figure  3  (after  Grossman  §  Melnik*)  shows  the  change  in  the 
boundary  layer  displacement  surface  when  a  propulsive  jet  replaces  a  sting. 
Figure  4  below  depicts  the  displacement  surface  location  for  two  different 
jet  velocities;  the  higher  the  jet  velocity  ratio  Uj/uQ,  the  more  rapidly 
the  boundary  layer  is  entrained  in  the  jet,  and  so  the  greater  the  suction 
on  the  afterbody.  Free  air  entrainment  by  the  jet  also  changes  the  static 
pressure  distribution  over  the  afterbody,  and  hence  the  forces  acting  upon 
it. 


Figure  4.  The  Boundary  Layer  Displacement  Surface  for  Two 
Different  Jet  Velocity  Ratios. 


This  jet-induced  drag  effect  was  noted  almost  as  soon  as  jet 

engines  were  applied  to  aircraft  propulsion.  Figure  5  shows  some  typical 

data  obtained  by  Riegels  and  Eggert  in  1944. 

3 

Kuchemann  and  Weber  observed  that  &Cp  the  change  in  the 

*  DC 

pressure  coefficient  Cp  varied  as  (uj/uq  -1)  and  was  practically  constant 
over  the  boattail,  so  that  the  pressure  and  drag  increments  due  to  entrain¬ 
ment  could  be  written  as 

ACp  =  0.01  (u./u  -1)  (1) 

BE  J  0 

0.01  [Uj/uQ  - 1  ]  [1  -  (r.j/rm)2] 

0.005 [P0/Pj ] [uo/Uj ] [ (rm/rj) 2  -  1] 
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Kuchemann  and  Weber  suggested  that  the  drag  increase  would  be 
typically  less  than  1%  of  the  thrust,  so  that  the  entrainment  effect  would 
be  unimportant.  This  is  no  longer  true  with  many  modem  aircraft  and 
missile  configurations,  of  course. 

There  are  two  alternative  methods  of  computing  the  change  in 

boattail  drag  due  to  entrainment  (Cnoc)*  The  most  time  consuming  is  to 

pfc 

compute,  iteratively,  the  pressure  distribution  over  the  boattail,  includ¬ 
ing  the  entrainment  effect  of  the  boundary  layer.  Without  entrainment, 
such  calculations  have  been  carried  out  by  many  investigators.  Recent 
examples  are  the  papers  of  Grossman  and  Melnik*,  Moulden,  Wu  and  Spring^, 
and  Chow,  Bober  and  Anderson'’.  But  such  studies,  while  they  are  impor¬ 
tant  to  gaining  a  better  understanding  of  the  underlying  physical 
processes,  are  not  at  all  suitable  for  correlating  experimental  data. 

And,  of  course,  additional  work  is  needed  to  introduce  a  satisfactory 
model  of  jet  entrainment. 

The  other  approach  is  to  employ  the  principle  of  linear  super¬ 
positioning.  If  it  is  possible  to  describe  the  flow  over  the  boattail  by 
a  suitable  distribution  of  singularities,  then  the  change  in  this  flow  due 
to  other  singularities  can  be  calculated  without  knowing  the  details  of 
the  original  ones,  or  the  details  of  the  corresponding  flow  field.  The 
limitations  of  this  approach  are  obvious.  The  displacement  surface  of 

the  boundary  layer  can  change  in  a  nonlinear  way  when  the  external  flow 
field  is  changed,  and  not  all  axisymmetric  bodies  can  be  described  by 
singularities  on  their  axis.  So  we  must  hope  that  the  errors  so  intro¬ 
duced  are  not  too  large  or  that  they  are  systematic  and  identifiable. 

It  seems  reasonable  to  suppose  that  jet  entrainment  can  be 
represented  by  a  sink  distribution  in  the  jet,  ideally  a  cylindrical  sink 
surface  located  at  the  jet  boundary.  But  since  we  do  not  know  the  entrain¬ 
ment  details  very  accurately,  a  line  sink  distribution  on  the  centerline 
axis  is  probably  a  reasonable  approximation.  But  before  attempting  to 
define  values  for  such  a  distribution,  we  need  to  describe  a  jet  in  a 
little  more  detail. 
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If  uQ  and  ikq  are  the  free-stream  and  jet  ("potential  core") 
velocities,  v  the  jet  radial  velocity,  and  x^  and  r  are  cylindrical  coordinates 
in  the  jet,  then  the  equations  of  motion  for  the  jet  are: 


3u. 


3u . 


u.  T— +  v  „ 
j  3Xj  3r 


=  1  —  (r 

r  3r  V  3r  / 


(3) 


!!i  +  |i  +  v  =  0 

3x .  3r  r 
3 


It  is  usual  to  express  the  virtual  kinematic  viscosity  e  in  the 
form  given  by  Prandtl's  mixing  length  theory* 

K 


e  =  r.u.  (tn  c  -  r.)|u.  -  u  I 

]  jo  2  0.5  1  jo  o' 


(4) 


where  Ky  is  the  (experimentally)  determined  velocity  mixing  coefficient 

rQ  5  is  the  radial  distance  at  which  the  local  velocity  is  the 
mean  of  the  core  and  free-stream  velocities 
[i.e.,  uQ  5  =  (u.o  -  uo)/2] 

r j  is  the  radius  of  the  "potential  core",  the  inner  edge  of  the 
mixing  layer. 

Typically,  K  is  correlated  against  the  Mach  number  correspon- 
*  ^ 

ding  to  Ug  j..  According  to  Smoot 


Ky  =  .047  for  0  <  MQ  5  <  0.2 


=  .028  for  0.7  <  MQ  5  <  1.3 

So  there  should  be  no  major  difference  between  low-speed  data 
and  transonic  entrainment.  Then,  following  Albertson,  et  al.  and  Kuchemann 
and  Weber,  the  entrainment  ratio  n  is  given  by 


*  The  Taylor  mixing  theory  gives  the  same  form;  Von  Karman's  somewhat  dif¬ 
ferent.  The  differences  are  not  important  here. 
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free  air  entrained  •>  /-ei 

n  ~  initial  jet  volume  flow  ^Xj  j?  o  jo 

(k  =  0.0415  at  low  Mach  numbers,  according  to  Albertson  et 

This  relationship  is  applicable  to  that  part  of  the  jet  close 
to  the  nozzle,  the  so-called  "zone  of  establishment"  where  it  still  has  a 
"potential  core"  of  velocity  UjQ.  This  limitation  will  be  justified  below. 

It  follows  from  equation  (5)  that  the  corresponding  axial  sink 
line  of  strength  q(xj)  per  unit  length  is 

q(Xj)  =  -Trkrj  (Uj0  -  uQ)  (ft3/sec.ft. ) 

the  minus  sign  reminding  us  that  it  is  a  sink  rather  than  a  source  line. 

From  Appendix  A,  the  velocities  induced  at  a  point  (x,r)  by  a 
constant  strength  sink  in  the  region  0  <  Xj  <  00  are 


(x  -  xx)  dx1 
[(x  -  Xj)2  +  (Br) 


213/2 


Br2  dXj 

7  rr"~7~Fj372 

(x  -  Xj)  +  (Br) 


Figure  6  shows  that  most  of  the  effect  on  Au  is  due  to  the 
first  five  diameters  or  so  of  the  jet,  thus  justifying  our  neglect  of  the 
far- jet  entrainment  functions.  [The  "zone  of  establishment"  typically  has 
a  length  of  roughly  9t^/\,  where  A  =  (Uj/uQ  -  1)].  With  q(Xj)  a  constant, 
equations  (6)  and  (7)  may  be  immediately  integrated  to  give 


q(xj 


Au  = 


/  5  2 

'(x  -  xx)2  +  (Brj 


42  +  (Br) 2 


q(xx) 


n  2 

x  +  (Br) Z 


as  Xj  +  “ 


Av  =  — 


q(Xj) 

4tt  r 


(x  -  Xj) 


7(x  -  xx)2  +  (6r)2  v/x2  +  ($r)‘ 


q(xx) 

4irr 


+  1 


L  / 


2  2 
'x*  +  (Br)z 


as  Xj^  -*• 


(9) 


Making  the  substitution  for  q(Xj): 


Au  = 


kr.  (u.  -  u  ) 

_J-  J° _ 9-1 

Jx  +  (6r)2 


(10) 


Av  = 


(uj0  ■  v* 


Jx2  ♦  (Sr)2 


+  1 


(ID 


These  are  plotted  in  Figure  7.  It's  clear  that  as  x^  •+• 


Au  -*■  0.  As  x^  -*■  ■»,  Au  •+  0  but 


Av  (k/2)(rj-/r)(Uj0  -  uQ) 
which  is  simply  the  flow  into  a  two-dimensional  sink,  of  strength  irkr - (u . 


uq) .  It's  also  clear  since  k  =  0(10  1)  that 


J  jo 


Au 


ft)’ 

ftf 


Thus,  from  equation  A16  of  Appendix  A: 


AC, 


Au 


BE 


-Me-1) 


-y 


'-2  „2*2 
x  *  S  r 


/ 


2  „2'2 

x  +  6  r 


(12) 


12 


/ 


Entrainment  Sink  Model. 


where  x  =  x/r^  r  =  r/r^ 

»  -  Ck/2Ku.0/u0  -  1)  .  U 

A  simple  numerical  example  of  equation  (12)  is  given  in  Figure 

8  for  a  cylindrical  afterbody,  r  =  r..  The  value  of  Cn  will  be  smaller 

m  j  P 

when  r  >  r. ,  of  course.  ^ 

m  j 

It's  now  possible  to  integrate  equation  (12)  over  an  afterbody 

in  order  to  determine  Cn  . 

U8E 

One  convenient  representation  for  an  afterbody  is  a  polynomial 

r  =  Tj  (1  +  Az  +  Bz2  +  .  .  . )  (13) 

where  z  =  -x-,  the  longitudinal  ordinate  measured  forward  from  the  jet  exit 
plane.  An  element  of  surface  is  2irr  ds  and  the  projected  area  2irr(dr/dz)dz. 
Thus,  the  drag  change  is 


d(ADE) 


1  2  _  dr  , 

2*>uo  ACPflP  2lTr  dldz 

pc 


dz 


(agd  ) 

UPE 


2 


r 

m 


r 


dr 

dz 


Substituting  for  ACp  from  equation  (12) 


d_(AC 

dz 


D 


2yr  dr/dz 

2  ri  ~n 

TmJz  *  B  r 


V 


L 


z 


m 


2 

(1  »  Az  +  Bz  +  ...) (A  +  2Bz  +  ...)  dz 
Jz  +  82  (1  +  Az  +  Bz2  +  ..)2 


(14) 
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which  can  be  numerically  integrated  for  any  number  of  terms.  For  the 
simplest  case  of  a  conical  afterbody: 


(1  +  Az)dz _ 

/  2  5  Y 

(1  +  0  A  )z  +  2A6  z  +  6 


where  A  =  tan  0,  the  tangent  of  the  cone  half  angle. 
This  is  conveniently  written  in  terms  of 


=  8A,  z  * 


r  -  1  dz 
A  »  df 


Then  6AC 


DBE  2*  / 

"  '  3  J. 


m  1  v  (1  4  -  2t  *  1 


'  3  I  a  *  *2> 


i  2  a  2  a 

(1  +  )r  -  2r  +  1 
m  m 


i  r  i  /  «  *  ♦2)i.  -  i  \  ,  i 

*  (J\  ,,2)3/2  lsinh  (, - * - /-si"h  »>J 


This  equation  is  plotted  in  Figure  9,  which  gives  AC^  for  any 
combination  of  rm,  6  and  M.  Some  specific  values  are  given  in  Figure  10. 
Note  that  the  "shape  factor"  «  =  ACp  /u  is  a  function  only  of  afterbody 
shape  and  freestream  Mach  number.  T^ioretical  values  of  ACp  for  a  parti¬ 
cular  conical  afterbody  are  plotted  in  Figure  11  for  three  different 
NPR's.  It  seems  clear  that  the  entrainment  drag  coefficient  is  less 
important  at  transonic  speeds. 

Using  the  approach  described,  it's  obviously  possible  to 

express  ACn  for  any  afterbody  as 
0E 

aCD  =  Ck/mu-o/Uo  -  1)*  (17) 
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tEAM  MACH  NUMBER 


t  for  Some  Representative  Conical  Afterbody 
tion  of  Mach  Number. 


j 

5 


\ 

j 

i 


Figure  11.  Theoretical  Entrainment  Drag  Increment  for  a  Conical  Nozzle 
Having  a  15°  Half  Angle.  r  =  2.0,  k  =  .0415 


where  the  shape  factor  4  is  given  by 


4  = 


r  m  , 

2  I  r(dr/dz)dz  2  / 

*2  r-= - = - 7  “  -2  I 

r  J  /  2  Q2  -2  r  . 

m  •'o  vz  +  6  r  ml 


Lm 


r  dr 


n  1^2 

sf  z  +  6  r 


(18) 


2  2-2 
Note  that  when  z  <<  $  r 


Br 


(r  -  1) 


(19) 


m 


Equation  (17)  now  explains  some  previously  puzzling  phenomena 
and  enables  us  (in  principle)  to  better  reduce  experimental  data  to  a 
common  denominator.  Figure  12  shows  some  of  the  Reubush^  data  plotted 
against  (Mj/Mq  -  1),  which  we  have  assumed  to  be  the  same  as  (UjQ/uo  -  1), 
since  the  jet  was  "cold."  The  raw  data  of  reference  8  shows  that,  for  a 
given  NPR,  reduces  with  increasing  Mach  number;  an  anomalous  result. 

But  when  plotted  as  shown  in  Figure  12,  we  see  that  this  is  explained  by 
the  reduction  of  ACL  with  increasing  Mach  number,  and  that  C  ,  the  value 

bE  Uftn 

for  zero  entrainment,  (given  by  the  intercepts  in  Figure  12)  increases  with 
Mach  number  as  we  would  expect. 

Figure  13  gives  a  more  detailed  plot  of  Cn  ,  for  another 
circular  arc  nozzle,  which  was  obtained  in  the  same  way.  Most  unfortunately, 
there  is  very  little  data  available  in  the  literature  for  subsonic  jet 
velocities,  so  that  it  is  not  possible  to  produce  parametric  curves  of 
Cp  ,  as  has  been  done  for  the  jet-off  case.  CD  for  some  NASA  circular 
ar^°boattails  is  plotted  in  Figure  14,  as  a  funi?ion  of  their  equivalent 
cone  angle,  and  considering  the  accuracy  with  which  the  raw  data  can  be 
obtained  from  the  original  references  and  the  extreme  paucity  of  subsonic 
jet  data  points  (usually  two;  never  more  than  three),  the  scatter  is  consi¬ 
dered  to  be  reasonable.  The  basic  data  for  these  boattails  is  given  in 
Table  1.  The  slope  of  the  data  in  Figure  12  gives  the  value  of  the  en¬ 
trainment  function  k,  since,  from  equation  (17) 

2ACL 


k  = 


BE 


4(u.  /u -  1) 
jo'  o 


(20) 
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BOATTAIL  DRAG  COEFFICIENT  CD 


FREE-STKEAM  MACH  NUMBER  M0 

Figure  13.  Ze  -o  Entrainment  CD  for  a  NASA  Circular  Arc  Afterbody  (NASA  TN 
P-7192),  DP/D..  =  O.S?  t/DM  *  0.8. 


TABLE 

:  l. 

DATA 

FOR  NASA  CIRCULAR 

ARC  NOZZLES,  FROM 

REFERENCES 

8  AND  11. 

°E 

£ 

c 

(Hi  -  ') 

dm 

M 

$ 

Deo 

Vuo  / 

k 

9f 

e 

c 

0.6 

1.0 

0.4 

0.404 

.0147 

.0114 

.0564 

22.07° 

11.03° 

0.6 

0.442 

.0184 

.0091 

.0412 

0.7 

0.475 

.02 

.0081 

.0341 

0.6 

1.5 

0.4 

0.346 

.0145 

.0084 

.0486 

14.81° 

7.41° 

0.6 

0.373 

.0171 

.0044 

.0236 

0.7 

0.396 

.0186 

.0019 

.0096 

0.7 

1.0 

0.4 

0.354 

.0107 

.0099 

.0559 

16.50° 

8.25° 

0.6 

0.388 

.0153 

.0048 

.0247 

0.7 

0.417 

.0167 

.0018 

.0086 

0.7 

1.5 

0.4 

0.305 

.0086 

.0073 

.0479 

11.04° 

5.52° 

0.6 

0.330 

.012 

.0044 

.0267 

0.7 

0.352 

.014 

- 

- 

0.7 

2.0 

0.4 

0.269 

.0086 

.0054 

.0401 

8.29° 

4.15“ 

0.6 

0.288 

.0077 

.0027 

.0188 

0.7 

0.305 

.0084 

.0014 

.0092 

0.5 

0.8 

0.4 

0.45 

.0243 

.0167 

.0742 

34.04° 

17.02° 

0.6 

0.496 

.0296 

.0154 

.0621 

0.8 

0.599 

.0382 

.0168 

.0561 

0.5 

1.0 

0.4 

0.419 

.0137 

.013 

.0621 

27.52° 

13.76° 

0.6 

0.458 

.0175 

.0103 

.0450 

0.5 

1.768 

0.4 

0.328 

.0068 

.0087 

.0530 

15.77° 

7.89° 

0.6 

0.351 

.018 

.0057 

.0325 

Average  k  =  .0548  for  M  =  0.4 

=  .0343  for  M  =  0.6 

=  .0154  for  M  =  0.7 
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The  shape  factor  <t> ,  obtained  by  numerical  integration  of  equa¬ 
tion  (18),  is  given  in  Figure  15,  and  the  derived  values  of  k  in  Figure  16. 
It's  encouraging  to  note  that  k  is  of  the  right  order  and  that  it  falls 
with  increasing  Mach  number.  The  increase  of  k  with  the  equivalent  cone 
half-angle  is  presumably  due  to  the  nonlinear  boundary  effects  mentioned 
at  the  beginning  of  this  discussion.  A  more  refined  analysis  would  recog¬ 
nize  that  u  <  uq  in  the  boundary  layer  leads  to  higher  entrainment  rates 
and  that  the  thicker  the  boundary  layer  (i.e.,  the  greater  the  equivalent 
half  angle),  the  more  important  this  effect  would  be.  A  more  refined  analysis 
would  also  compute  the  important  change  in  the  boundary  layer  displacement 
thickness  with  increasing  entrainment. 

NSRDC  boattail  10  has  no  base  area  and  so  can  be  compared 
with  the  theory  as  shown  in  Figure  17.  The  value  of  k  =  0.02  was  extrapolated 
from  Figure  16,  which  clearly  needs  to  be  expanded.  The  value  of  4>  was  read 
from  Figure  9,  even  though  boattail  10  is  not  a  simple  cone,  because  the 
fairing  into  the  body  is  furthest  removed  from  the  jet,  and  the  error  should 
therefore  be  tolerable. 


PLUME  INTERFERENCE 

As  indicated  in  Figure  18,  plume  interference  is  characterized 
by  an  increase  in  static  pressure  on  the  boattail.  If  the  jet  diameter  is 
smaller  than  the  base,  the  length  of  the  separation  streamline  tends  to 
partially  insulate  the  boattail.  Thus,  when  there  is  a  substantial  base 
area  "insulating"  the  boattail  from  the  jet,  the  plume  effect  on  boattail 
drag  can  be  negligible  for  moderate  values  of  NPR. 

The  first  step  in  understanding  the  phenomena  is  to  estimate 
the  shape  of  the  plume.  To  do  this,  we  first  determine  the  initial  turning 
angle  ,  defined  in  Figure  19. 

It's  been  observed  by  many  workers  that  the  initial  jet  turning 
angle  6^  is  not  much  affected  by  the  free  stream  flow  conditions,  so  that 

6.  =  v.  -  v  +  a.,  (21) 

j  INN 
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o  &2  0.4  0.6  O.e  1.0 

FREE  -  STREAM  MACH  NUMBER 

Figure  15.  The  Shape  Factor  $  for  Some  NASA  Circular  Arc  Boattails. 


Figure  16.  Variation  of  Entrainment  Ratio  k  with  Free-StTeam  Mach  Number 
and  Equivalent  Cone  Angle  for  the  NASA  Circular  Arc  Nozzles. 
(The  precise  form  of  the  variation  is  not  known  at  the  present 
time. ) 
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Figure  17.  A  Comparison  Between  Entrainment  Theory  and  DTNSRDC 
Measurements  on  Body  70010  at  M  ■  0.7. 


Figure  18.  The  Magnitude  of  Plume  Interference  on  the  Boattail  Depends 
on  the  Ratio  of  Jet  to  Base  Diameters. 
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Figure  19.  Initial  Jet  Turning  Angle 


where 


where  Mg 


Uj  =  the  Prandtl-Meyer  expansion  angle  for 
the  plume  boundary  Mach  number  Mn 

D 


(mbz  -  1) 


-tan 


-1 


the  Prandtl-Meyer 
the  nozzle  exit  Mach 


expansion  angle  for 
number  M^ 


K2/*)  - tan‘1  A 


and  is  specified  by  the  nozzle  design 


If  the  flow  were  one -dimensional  the  maximum  cross-sectional 
area  (Am)  of  the  plume  would  be  given  by 


\id  ■  £  [rfr  (>  *  ^  «b2)  ] 


v  t  1 
2(Y  -  1) 


Y  ♦  1 
2  \  2(Y  -  1) 


y  ♦  l 

2Y 


W . (fr-f 

/W(£*)  V  - .] 


where  A  is  the  sonic  throat  area 

Equation  (25)  is  plotted  in  Figure  20  for  typical  pressure  ratios. 
It  turns  out  (as  shown  in  Figure  21)  that  for  zero  nozzle  angle  the  ac¬ 
tual  area  ratio  is  quite  precisely  the  square  of  this  quantity,  for  NPR's 
less  than  about  20.  A  simple  linear  approximation  is  also  shown  in  Figure  21. 

In  experimental  measurements  of  plume  boundaries,  the  static 
pressure  p£  at  the  jet  exit  is  often  employed.  This  is  related  -  again 
one  dimensionally  -  to  the  jet  exit  Mach  number  by 


(‘♦V-  »B2) 


Y 

Y  -  1 


Pe  .  Pe  fit 
p-  Pjt  p» 


£(>- 


1  \2) 


Y 

Y  -  1 


Underexpanded  jet  boundaries  have  been  studied  by  many  workers. 

Of  particular  note  are  the  works  by  Love  et  al.^  and  Vick  et  al.,^  both 
of  which  give  additional  important  references.  For  the  moderate  nozzle 


pressure  ratios  of  interest  here  it's  generally  concluded  that  the  jet 
profile  is  adequately  described  by  a  circular  arc  between  the  nozzle  exit 
and  the  maximum  diameter,  as  indicated  earlier  in  Figure  19.  Unfortunately 
the  equation  for  a  circular  arc  is  relatively  complicated  to  integrate 
when  using  slender  body  theory.  Accordingly,  we  here  employ  the  second 
order  polynominal 


r  =  r,  +  6 .x.  - 

3  3  3 


where  x .  2 (r  -  r . ] 

jm  =  m  jJ 

S. 

J 


For  initial  jet  expansion  angles  of  0  -  30° ,  this  gives  a 
profile  which  is  essentially  identical  with  a  circular  arc. 

From  the  theory  of  Appendix  A  we  can  now  compute  the  velocity 
pertubations  due  to  the  plume  and  the  resulting  pressures  on  the  boattail. 
From  equation  A. 18 


q(x.)  *  2ir  un  T 


o  dx. 


=  2ir 


Some  typical  q(x^)  variations  are  plotted  in  Figures  22  and  23. 
For  small  values  of  NPR 


q(Xj)  =  2-n  uq  r^  6^  (1  -  Xj) 


(xj  ■  W 


is  not  a  bad  approximation  to  equation  (29).  It  becomes  worse  as  NPR 
increases,  but  then,  so  does  the  slender  body  approximation  itself.  So 
although  our  discussion  is  initially  more  general,  we  shall  return  to 
this  approximation  for  an  engineering  approximation  to  the  plume  effect. 


4  7T  uft 


0.3 


Figure  22. 


The  Source  Distribution  Coefficient  for  Various  Nozzle 
Pressure  Ratios.  Sonic  Nozzle,  y  -  1.4 
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Figure  23.  The  Source  Distribution  Coefficient  for  Various 
Nozzle  Pressure  Ratios  as  a  Function  of  Nozzle 
Radius.  Sonic  Nozzle,  y  =  1.4 
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The  integration  is  most  conveniently  carried  out  if  we  use  the 
dummy  variable 


Then  r 


Z  *  X  -  X. 

3 

so 

X.  =  X  -  z 

3 

(Note  that  z 

differs  by 

Br  from  the  value 

=  r.  +  6.  (x 
3  3 

-  z)  -  6i 

2  2 
(xz  -  2xz  +  z^) 

*■>. 

■  b  *  v  (- 

jut J 

.  -  4; 


2x . 


jm 


A  -  Bz  -  Cz 


say 


(31) 


q(z)  =  2tt  uq  r  =  -  2tt  uQ  (A  -  Bz  -  Cz2)  (B  +  2Cz) 


when  a 


2  7 

=  - 2tt  uq  (a  ♦  bz  +  ez  +  fzJ)  say 
-  AB-  Vj.  [*j  *  (1  - 1  *>]  ”  - 


(32) 


(B2  -  2AC)  =  -6^ 2(1  -  x)2+  S 


3BC  =  -  |  (1  -  *) 

x. 

jm 


2  1 

2C  =  -  2  ^“2 

*jm 


where 


x/x 


jm 


ri  '  r)/Xjm 
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We  can  now  determine  the  velocity  pertubations  from  Appendix  A 
From  equation  A4,  noting  that  dz  =-dXj  the  axial  pertubation  is 

z. 


u 

u 


/  2  2 
laz  +  bz  +  ez  + 

P  *  .V) 


*‘4) 


dz 


1 _  (a  +  bz  -  eg2r2  -  f?62r2) 


£77~t2 


*  Vz^+gr  (e  +  jfz) 


(b  -  j  f$2r2)  log  (z  +  J  z2  *  g2r2) 


(33) 


From  equation  A5,  the  radial  velocity  pertubation  is  given  by 
•z. 


v  g2r  [  V  +  bz  4  ez2  4 

=  2  l  JJTsJW2 


fz3\  dz 


2_  r 


e_r 

2 


Lv£777 


(az 
62r2 


2  2 

b  -  ez  +  fg  r 


+  f  J2  +  g2r2  +  e  log  (  z  +  >/z2  4  g2r2  )  j 


(34) 


While  one  could  formally  proceed  to  expand  these  equations,  the 
form  given  is  most  suitable  for  numerical  evaluation.  The  limits  are 
Zj  =  x  -  x.  and  zq  *  x.  While  this  "exact'  solution  may  be  of  research 
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interest,  its  usefulness  is  limited  by  the  slender  body  approximations 
employed  to  define  q(Xj).  It  cannot  be  accurate,  or  even  close,  when  6^ 
is  large.  We  may  as  well  employ  equation  (30)  therefore,  so  that 

q(z)  =  2ir  u  r.  6.  (l  -  X  ~-Z  )  (35) 

o  i  i  \  x.  / 

J  J  jm  ' 

.'.a  =  -  r..  6^  (1  -  M)  (x  =  x/x^,  etc.) 


and  C 

P 


Again,  because  of  limited  accuracy,  we  may  as  well  neglect 

2 

(v/uQ)  in  relation  to  Au/uq.  Then  we  find  that  the  change  in  Cp  due  to 
the  plume  is 


(38) 
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But  this  introduces  the  additional  terms  in  6.  and  r  /r., 

1  m'  j 

which  depend  upon  p^+/pm  ,  so  for  a  simple  overview  of  the  plume  effect, 
equation  (36),  as  presented  in  contour  in  Figure  24,  is  often  preferable. 

It  is  clear  from  both  forms  of  the  equation  that  a  simple 
integral  (like  4>  in  the  section  on  entrainment)  cannot  be  derived  to 
account  for  the  plume  effect.  The  change  in  boattail  drag,  due  to  the 
plume,  will  be 


AC 


Dp 


r . 
1 


m 

AC 


2  r  dr 


where  AC 


6. 

3 


& 


r  "ft 

h  n2-2 

x  +  6  r 

log 

x  -  1  +  V(x  -  1)2+  B2?2  ( 

■i 

-  A 
"  / 

i  *  A2  ♦  b2;2  Jj 

(39) 


and  both  6.  and  rm  are  functions  of  P^/po,  ,  as  given  earlier. 

So  even  this  very  approximate  formulation  must  be  numerically 
integrated.  That  being  so,  one  might  as  well  employ  the  complete  slender 
body  equations  to  compute  ACpp.  This  has  been  done  in  a  Fortran 
conversational  program  PLUME3  listed  in  Appendix  B.  The  program  also 
includes  the  entrainment  effects  discussed  earlier,  equation  (17).  A 
comparison  of  the  program  output  with  experiment  is  .given  in  Figures  25 
and  26. 

As  might  be  expected,  there  is  some  disagreement  in  the  absolute 
values  of  C^,  particularly  close  to  the  jet,  because  our  inviscid  flow  theory 
neglects  the  boundary  layer  displacement  thickness.  But  the  change  in  C^ 
due  to  the  plume  effect  appears  to  be  correct. 
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Figure  24.  Contours  of  AC  in  x  and  r.  Ac  is  the  change  in 

JP.  PP 

?.  . 

J  J 

boattail  pressure  due  to  the  plume  shape. 

x.  *  downstream  distance  from  the  jet  exit  to  the  maximum 
•|m  plume  diameter. 

-x  =  upstream  distance  from  the  jet  exit  plane, 
r  =  local  boattail  radius. 

B  =  >/l  -  M2 

<5j  =  jet  turning  angle  at  the  nozzle. 
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Figure  25.  A  Comparison  Between  Theory  and  Experiment  for  M  ■  0.4 

Data  of  Reubush  and  Runkel.^  Configuration  3  afterbody. 
(Theory  from  PLUME  3  program.) 
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Figure  26.  A  Comparison  Between  Theory  and  Experiment  for  M  *  0.9. 

Data  of  Reubush  and  Runkel.11  Configuration  3  afterbody. 
(Theory  from  PLUMES  3  Program.) 


SUMMARY  AND  CONCLUSIONS 


A  theoretical  framework  has  been  presented  for  the  prediction 
of  the  change  in  boattail  drag  with  nozzle  pressure  ratio  (NPR)  and 
agreement  with  experiment  seems  good  in  those  few  cases  where  it  was 
possible  to  check.  The  most  severe  restrictions  are 

•  No  base  area 

•  Attached  flow  with  negligible  boundary  layer 
(BL)  thickness. 

The  second  of  these  would  be  fairly  simple  to  relax,  by  incor¬ 
porating  axisymmetric  BL  equations.  We  would  then  have  a  tool  which  per¬ 
mitted  the  optimum  shape  for  particular  internal  volume  constraints  to  be 
determined  on  the  computer.  Perhaps  more  importantly,  it  would  enable 
unsatisfactory  shapes  to  be  identified  and  modified  before  testing. 
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APPENDIX  A 

DERIVATION  OF  THE  SUBSONIC  SLENDER  BODY 
EQUATIONS  FOR  A  BODY  AT  ZERO  INCIDENCE 


(In  this  appendix,  the  perturbation  velocities 
Au,  Av  are  written  as  u,  v  for  conciseness. 
Also,  tp  is  the  familiar  stream  function 
<f>  is  the  velocity  potential) 
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The  Representation  of  a  Body  by  Singularities 

A1 

As  was  first  suggested  by  Rankine  the  axi-symmetric  flow  about  large 
classes  of  bodies  can  be  represented  by  distributions  of  sources  and 
sinks  (negative  sources)  along  their  axes.  If  the  total  strength  (Q)  of 
the  sources  and  sinks  is  zero,  the  streamline  representing  the  body  will 
be  closed;  if  not,  it  will  be  semi-infinite,  the  extension  to  infinity 
being  downstream  if  the  total  strength  Q>0,  and  upstream  if  Q<0. 

The  velocity  potential  <j>  and  stream  function  for  a  single  source  of 
strength  Q  (ft3/sec)  at  the  origin  in  terms  of  cylindrical  coordinates 
(r,x)  is 


4irv/r2  +  x2 


and  the  stream  function  is 


♦ 


-Qx 


.  [~2  2 
4ir  y/r  +  x 


One  can  also  (or  alternatively)  employ  higher  order  singularities,  when 
convenient.  By  differentiating  the  above  equation  with  respect  to  x  we 
obtain  the  values  of  $  and  i Ji  for  a  source  doublet:  viz 


-Qx 


3/2 


4> 


4tt 


12 
>  r 


3/, 


+  x 


Further  differentiation  will  yield  singularities  of  higher  order  at  the 
origin. 

More  general  classes  of  body  shape  may  be  calculated  using  source  or 
vortex  rings  or  discs,  discrete  or  continuous,  either  "buried"  below  or 
at  the  stream  surface  which  represents  the  body.  In  general,  these  more 
sophisticated  singularities  involve  somewhat  more  complex  analysis 
(particularly  when  the  flow  is  compressible)  and  these  complications  are 
not  necessary  for  the  present  problem* 
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The  General  Distributed  Source  Body  in  Subsonic  Flow 

Small  pertubation  theory  in  subsonic  three-dimensional  flow  gives  the 
velocity  potential  of  a  point  source  located  on  the  x  axis  as  A2 


<f>  = 


m 

4tt 


1 

\/(x  -  XjJ2*  B2r2 


(Al) 


,  2  2  2 
where  r  =  y  +  z 


m  =  the  source  strength  (negative  for  a  sink) 

M  =  the  free  stream  Mach  number 

Xj  =  location  of  the  source  on  the  x  axis 


The  velocity  perturbations  u  and  v  (in  the  x  -  direction  and  radially 
normal  to  it)  associated  with  the  source  are  therefore 


u 


v 


9<(> 

3x 


3<t> 

3r 


(x  -  Xj) 

[(X  -  x x ) 2  +  S2r2J3/2 

32r _ 

(x  -  Xj)2  +  32r2j  3^2 


(A2) 


(A3) 


For  a  continuous  distribution  of  sources  q(x,)  between  x  =  a  and  x.  = 
it  follows  that 

m  =  q(x^)dxj 


u  = 


j_  f 

4tt  J 


q(xJ 


(X  -  Xj) 


[(X  -  xp2  +  e2r2 


3/2  dxl 


(A4) 


v  = 


/ 


2_  /  q(Xj)  r  -g  r - 

4lT  [(x  -  Xx)2  +  e2r2] 


3/2  dxl 


(A5) 
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Let  z 


-L_  ) 

Br 


(A6) 


(A7) 


Integration  by  parts  enables  us  to  express  the  velocity  components  as  a 

series  in  d\(z)  plus  a  residual  integral.  The  first  few  terms  are 
dzn 
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In  other  words,  the  flow  field  is  defined  by  the  nature  of  the  source 
line  at  each  end.*  If  the  third  derivative  of  the  source  strength  is 
zero  at  its  ends,  equations  (A8)  and  (A9)  give  the  flow  field  explicitly; 
and,  of  course,  the  equations  can  be  extended  to  higher  order  derivatives 
by  further  integration  by  parts. 

In  passing,  it's  of  interest  to  note  that  the  simplest  possible  body  of 
this  type  (apart  from  the  limit  case  of  a  sphere)  is  given  by 


*  =  %  '  <lix 


n 


(A10) 


for  the  body  to  be  closed 
-xb 


=  J qdXj  =  qQxb  -  \ 


x,  =  2q 
b  no 


f^b  =  0 

qo  -  *l*b  =  -% 


From  equations  (A8)  and  (A9) : 


u  = 


4tt 


l^x2  ♦  g2r2  •  / 


v  = 


4Ttr 


/  2 

a2  2 

V  (X  -  xb) 

+  B  r 

2  10CT 

x  +  /c2  ♦ 

D 

jx  '  xb)  + 

(x  -  xb) 

(All) 


_  Bx  +  +  2_ 

\/x2  +  (lV  \/(x  -  xb2)2  +  g2r2  Xb 


2  J~2  ~2  2 

+  —  vx  +  g  r 


-  b  /(*  -  xb) 


2  a2  2 
+  g  r 


(A12) 


*This  is  the  reason  why  not  all  body  shapes  can  be  represented  by  this 
type  of  singularity. 
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v  =  o  and 


At  the  midpoint,  x  =  jxb. 


~~  log 
2™b 


4  +  (28r/xb)2  +  1 
4  +  (28r/xb)2  -  1 


(A13) 


For  any  body,  the  boundary  condition  at  the  surface  is 


v  dr  _  dr  dz 

u  +  u  dx  ~  dz  dx. 


1  dr 
8r  dz 


(A14) 


dr 

So  since  is  known,  we  can  in  principle  solve  for  q(z),  the  source 

distribution  required  to  give  r(z).  In  practice  it  is  easier  to  seek 
the  desired  shape  iteratively. 

The  resultant  velocity  is 


V  =  V(u  +  u)2  +  v2 
o 


(A15) 


Thus,  from  Bernoulli,  the  local  static  pressure  p  is  given  by 

12  1  T  2  2T 

Poo  +  2^00%  =  P  +  ^-L^O  +  U)  +  V  J 

p  -  P  r  2  2  1 

■■■v  c?  '  -  -iy  ] 


(A16) 


/  u  \^* 

Note  that  the  (—  J  term  can  generally  be  disregarded  in  evaluating  (A16), 
uo 

but  not  always  /  v_\2,  This  is  because  / v_\  is  often  much  larger  than  /  u  \ 
'un/  'uJ  'u_/ 


The  Slender  Body  Approximation 

"Slender  body"  theory  assumes  that  the  body  is  so  slender  that 
Br<<  (x  -  Xj)  and  u<<uQ.  Thus,  from  equation  (A14) 


dr  „  v 
dx  "  ur 


(A17) 
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Following  the  early  work  of  investigators  such  as  LaitoneA4,  WardA^,  NeumarkAC> 
and  Adams  and  SearsA? 


dS 

dxi 


q(xl)  “  u0  ^  =  2nurtr 


dr 
ox  dxj 


(A18) 


where  S  is  the  local  body  cross-sectional  area  irr2 ,  and 
assuming  that  there  are  no  local  discontinuities  in  dr/dx^,  then 


q(z)  =  .  2,mO/dr 

8  'dz 


(i)  (“  - 1) 


(A19) 


Using  (A18)  in  conjuction  with  (A4)  and  (A5)  or  (A19)  with  (A6)  and  (A7) 
solutions  for  u  and  v  are  readily  obtained.  Analytic  solutions  are 
possible  for  simple  body  shapes,  but  numerical  integration  of  (A4)  and 
(A5)  is  extremely  simple,  and  is  much  to  be  preferred  for  general  purposes. 


The  Tapered  Sting  Problem 


A3 


This  was  first  essayed  by  Tunnell  using  the  slender  body  equations  of 
Laitone.AV  These  equations  neglect  the  radial  velocity  pertubation  (v) , 
which  we  will  include. 


For  o<x<t,  r*r+x  tan  9.  To  calculate  the  effect  of  the  sting 
taper  on  the  model  pressure  distribution,  we  need  to  solve  for  C  in 
the  region  x  <  o. 


This  is  directly  related  to  a  more  general  problem;  that  of  computing  the 
inviscid  pressure  distribution  over  the  surface  of  a  conical  afterbody, 
and  the  change  in  pressure  it  induces  elsewhere.  In  the  notation  of  the 
sketch  above,  the  solution  is 


u 

u 

o 


rQ  (rQ  ♦  l  tan  9) 

+  $2r2  /(x-Jt)2  +  B2r2 


-  tan  0  log 


(A20) 


2  tan  0 


2 

x-H  +  B  r  tan  0 
/(x-t)2  +  B2r2 


x  ♦  B2r2  tan  0  ) 

/x2  7  B2r2  j 


(A21) 


Cp  =  -  [2 (u/uQ)  +  (u/uq) 2  +  (v/uq)2] 


(A16) 


(When  using  these  equations  for  a  conical  boattail,  it  should  be  remem¬ 
bered  that  0  will  be  negative,  in  contrast  to  the  usual  convention.)  If 
we  write  r  =  r  ,  equation  A20  is  identical  with  Tunnell's,  as  would  be 
expected. 

Tunnel 1  obtained  a  solution  for  r  =  rQ,  since  he  was  only  interested  in 
the  effect  of  the  sting  taper  on  base  pressure.  Figure  A1  gives  a  com¬ 
parison  of  his  equation  with  those  above,  and  his  experimental  data.  Both 
analyses  neglect  the  boundary  layer  displacement  thickness  over  the  model 
and  the  sting.  In  particular,  the  boundary  layer  will  reduce  the  effec¬ 
tive  angle  0  of  the  sting  taper,  and  increase  rQ;  both  of  which  would  re¬ 
duce  the  thebretfcdl  "values  of  SC  slightly,  especially  in  the  vicinity 
of  the  discontinuity.  ^b 

Since  the  sting  taper  was  always  present  in  Tunnell's  experiments,  we  have 
no  means  of  knowing  the  "baseline"  value  to  subtract  in  order  to  get  AC  , 
the  change  due  to  the  taper.  But  the  theory  should  be  most  accurate  ^ 
when  the  taper  is  furthest  removed  from  the  base;  so  this  value  of  AC  is 
used  to  compute  the  baseline  value. 

Note  that  by  integrating  equation  A16  over  a  boattail,  we  can  determine 
the  effect  of  a  tapered  sting  on  Cp  as  well  as  Cp  . 

Figure  A2  compares  equations  A20-A21  with  measured  pressure  distributions 
over  a  conical  boattail.  Agreement  is  rather  poor  in  this  case.  Near  x  =  0 
the  discontinuity  in  r  causes  a  large  suction  "spike"  which  slender  body 
theory  cannot  predict.  (Shape  discontinuities  cannot  be  produced  by  any 
combination  of  singularities  on  the  axis  of  symmetry.)  And  the  combined 
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effects  of  jet  entrainment  and  boundary  layer  growth  (mostly  the  latter) 
prevent  the  full  recompression  being  achieved  near  the  exhaust  plane. 


Comparison  Between  Slender  Body  Theory  and  Exact  Solutions 

Several  authors  have  compared  slender  body  theory  with  Lamb's^  exact  solu¬ 
tion  for  a  spheroid.  The  horizontal  velocity  disturbance  at  the  point  of 
maximum  thickness ,  for  example, is  given  by  Lamb  as 


where 


u 

u 

o 


max 


<52  (log  e  +  2f) 
2f  -  62  log  e 


f 


e 


7l  -  <Va>2 

1  +  f 
1  -  f 


(a  =  semi-major  axis) 


Slender  body  theory,  on  the  other  hand,  gives 


u~  -  T  (Va>2  l0* 

o 

max 

The  corresponding  values  of  u/uq  are  tabulated  below. 


rm/a 

Exact 

solution 

Slender  Body  Theory 

.0001 

8.9  x 

ID'8 

9.9  x 

10'8 

.001 

6.6  x 

10"6 

7.6  x 

lO"6 

.01 

4.3  x 

10~4 

5.3  x 

10'4 

.1 

.0207 

.03 

.2 

.0591 

.0925 

1  +  A  *  (rm/a)2 
-  1  +  A *  (r  /a)2 


The  discrepancy  here  is  due  to  the  fact  that,  no  matter  how  low  its  thick¬ 
ness  to  chord  ratio,  the  ends  of  a  spheroid  are  not  "slender"  and  introduce 
large  errors  into  the  equation.  So  the  comparison  is  not  really 
valid.  A  more  appropriate  comparison  would  require  the  exact  solution  to 
be  for  a  body  with  pointed  ends,  or  for  a  nondiscountinuous  change  in  the 
diameter  of  an  infinite  body. 
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APPENDIX  B.  PROGRAM  "PLUME3" 


This  program  uses  equations  (17),  (33)  and  (34)  to  compute  the 
pressure  distribution  over  circular  arc  boattails.  The  arc  size  can  be 
anywhere  between  infinity  (which  gives  a  conical  boattail,  of  course) 
and  the  minimum  which  fits  tangentially  to  the  forebody  and  gives  the 
correct  surface  angle  8  at  the  nozzle  exit  plane.  In  retrospect,  it 
would  have  been  better  to  have  used  a  power  series  representation,  since 
slender  body  theory  breaks  down  when  the  surface  is  discontinuous.  Thus 
only  the  tangential  arc  case  is  really  usable  in  this  program.  However, 
the  modifications  required  to  generalize  the  boattail  shape  are  relatively 
st rai ght - forward . 

The  program  is  "conversational,"  a  typical  output  being  shown 
in  Table  B.l.  The  only  difference  between  these  two  runs  is  that  k  =  0 
in  the  first  and  k  =  .05  in  the  second.  The  inputs  are 
R-M  =  r  ,  the  maximum  body  radius 
R-J  =  r . ,  the  initial  jet  radius 
L-BETA  =  the  length  of  the  afterbody 

P  =  NPR,  the  nozzle  pressure  ratio 
GAMMA  =  y»  the  ratio  of  specific  heats  in  the  jet 
K  =  k,  the  entrainment  coefficient 
VS  =  a,  the  speed  of  sound  in  the  jet 
M  =  M,  the  free  stream  Mach  number. 

Table  B.2  lists  the  program. 


Table  B.l. 


DO  YOU  WANT  A  NEW  RUN?  TYPE  YES  OR  NO. 
•YES 

NEW  VARIABLE?  !M».9 
NEW  VARIABLE?  !NOrl 
INPUT  VARIABLES  FOR  THIS  RUN 


R-M 

R-J  L-BETA 

P 

GAMMA 

1 

.5  3.536 

2.02 

1.4 

K 

VS  M 

N 

N1 

0 

1117  .9 

40 

200 

KSI 

1 

CALCULATED 

R-P 

PARAMETERS  FOR  THIS  RUN 

THETA-J  BODY  RADIUS 

PLUME  RADIUS 

PLUME  TERM 

.501 

212  . 970882E-02  12.7532 

10618.5 

106.627 

X 

U/U-0 

v/u-o 

C-P 

C-D  INT 

0 

. 620742E-01 

- .  142390E-01 

-.128204 

- . 313574E-07 

.3536 

.800171E-01 

- . 321 1 10E-01 

-.167468 

. 194095E-02 

.7072 

.870689E-01 

- . 554903E-01 

-.184798 

•  76221SE-02 

1.0608 

.850697E-01 

- . 807544E-01 

-.183897 

.  168544E-01 

1.4144 

• 759378E-01 

-.10685 

-.169059 

• 285567E-01 

1.768 

. 601045E-01 

-.133474 

-.141637 

. 410201E-01 

2.1216 

. 368033E-01 

-.160467 

-.100711 

. 520145E-01 

2.4752 

• 353787E-02 

-.187509 

- . 422479E-01 

. 587809E-01 

2.8288 

- . 467444E-01 

-.213283 

. 458143E-01 

. 577164E-01 

3.1824 

-.13956 

-.229367 

.207034 

• 427128E-01 

3.536 

-.314962 

-.144381 

.509877 

. 131749E-02 

64 


rj  ro 


DO  YOU  WANT  A  NEW  RUN?  TYPE  YES  OR  NO. 
!  YES 

NEW  VARIABLE?  IK..05 

NEW  VARIABLE?  !NO.I 

INPUT  VARIABLES  POR  THIS  RUN 


R-M 

R-J 

L-BETA 

P 

1 

.5 

3.53A 

2.02 

K 

VS 

M 

N 

500000E-01 

1117 

.9 

40 

KSI 

1 

CALCULATED  PARAMETERS  FOR  THIS  RUN 


R-P 

THETA-J 

BODY  RADIUS 

PLUME  RADIUS 

.501212 

. 970882E  02 

12.7532 

1061 8.5 

X 

U/U-0 

V/U-0 

C-P 

0 

. A29333E-01 

-. 14A134E-01 

-.130041 

.  353A 

.80A8A3E-01 

- . 327250E-01 

-.168954 

.7072 

. 875221E-01 

- . 5A1992E-01 

- . 1B58A3 

1 .0A08 

. 853983E-01 

- . B15108E-01 

-.184733 

1.4144 

.7A1925E-01 

-. 107A48 

-.169778 

1 .7A8 

• A031 15E-01 

-.134319 

-.142302 

2.121A 

. 3A9774E-01 

-.1A1372 

-.101363 

2.4752 

. 3A8782E-02 

-.188494 

- • 429191E-01 

2.8288 

-  •  4AA12BE-01 

-.214377 

. 450951E-01 

3.1824 

-.139443 

- . 230A2 

.206256 

3.53A 

-. 314S5A 

-.145877 

.509298 

DO  YOU  WANT  A  NEW  RUN?  TYPE  YES  OR  NO. 
!  NO 

2290  EXIT 
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GAMMA 

1.4 

N1 

200 


PLUME  TERM 
106.627 


C-D  INT 

-.3180A5E-07 
. 19A058E-0 
, 767915E-0 
. 1695A4E-01 
. 287091E-01 
. 412281E-01 
.522841E-01 
. 591 189E-01 
• 581303E-01 
. 432097E-01 
. 1888A3E-02 


/ 


W  M 


PAGt 


1 


DI ALUuM  PLUME3 


ID#380  06/28/77 


Table  B.2. 

10  DIh  A*( 15) 

20  READ  R1  .R2.L.P.G9.K5.M.N.N1  .K9.A0 
*’  30  DATA  1.. 5. 3. 536. 2. 02.1.4.. 05. .4.40.200. 1.1117 
40  FOR  1=1  TO  13 

50  READ  A*(I> 

60  NEXT  I 

70  DATA  R-M.R-J.L-BETA.P.GAMMA.K.VS.M.N.Nl »KSI .YES. NO 
80  PRINT  ’DO  YOU  WANT  INSTRUCTIONS?  (TYPE  YES  OR  NO)'? 

90  INPUT  B* 

100  IF  Bi=A$(13>  GO  TO  600 
110  IF  B*  =  At (12)  GO  TO  130 
120  GO  TO  80 

130  PRINT  'DEFAULT  VALUES  FOR  ALL  VARIABLES  HAVE  BEEN  ASSIGNED.' 

140  PRINT  'THEY  WILL  BE  LISTED  AFTER  THESE  INSTRUCTIONS.' 

150  PRINT  'THEN  YOU  WILL  BE  ASKED  WHETHER  YOU  WISH  TO' 

160  PRINT  'CHANGE  A  VARIABLE.  AT  THAT  POINT  TYPE  IN' 

.170  PRINT  'THE  VARIABLE  SYMBOL  FOLLOWED  BY  A  COMMA  AND  THE' 

180  PRINT  'VALUE  YOU  WISH  ASSIGNED  TO  IT.  IF  YOU  TYPE' 

190  PRINT  '(NO,  ANY  NUMBER).  THE  RUN  WILL  BEGIN.  THE  VARIABLE' 

200  PRINT  'SYMBOLS  ARE  AS  FOLLOWS » ' 

210  PRINT 

220  PRINT  TAB(7)i ' R-M ' » TAB ( 20 ) » 'IS  THE  MAXIMUM  RADIUS  OF  THE  BODY' 

230  PRINT  TAB(7)»"R-J*»TAB(20)»*IS  THE  RADIUS  OF  THE  JET* 

240  PRINT  TAB(7) » 'L-BETA* .TAB  (20) »  *  IS  THE  HORIZONTAL  DISTANCE  ' 

250  PRINT  TAB(20) f 'BETWEEN  THE  POINTS  AT  WHICH  R-M  AND  R-J' 

260  PRINT  TAB (20) ? 'ARE  MEASURED.' 

270  PRINT  TAB(7)»'P'»TAB(20)»*IS  THE  PRESSURE  RATIO  AT  THE  JET' 

280  PRINT  T AB ( 7 ) » ' GAMMA 'fTAB(20)»*IS  THE  RATIO  OF  SPECIFIC  HEATS* 

290  PRINT  TAB( 20 ) ? 'AT  THE  JET* 

300  PRINT  TAB(7) »  *K' {TAB (20) > ' IS  THE  ENTRAINMENT  COEFFICIENT* 

310  PRINT  TAB(7)»'VS'»TAB(20)?*IS  THE  VELOCITY  OF  SOUND  AT  THE* 

320  PRINT  TAB(20>? 'STAGNATION  TEMPERATURE' 

330  PRINT  TAB(7).'M*»TAB(20)»*IS  THE  FREE  STREAM  MACH  NUMBER' 

340  PRINT  TAB ( 7) »  *N“ {TAB ( 20 ) . *  IS  THE  NUMBER  OF  STATION  POINTS  (X)' 

350  PRINT  TAB(20) i 'ALONG  L-BETA' 

360  PRINT  TAB (7)»"N1*»TAB(20)»'IS  THE  NUMBER  OF  INTEGRATION  POINTS' 

370  PRINT  TAB(20>  » "ALONG  L-BETA* 

380  PRINT  TAB( 7 >»*KSI*»TAB(20)?*I5  THE  FRACTIONAL  SAGITTAL  DISTANCE' 
400  PRINT  TAB(20) i 'BETWEEN  A  STRAIGHT  LINE  BODY  CURVE  AND  THE  DESIRED’ 
410  PRINT  TAB (20) ? 'CIRCULAR  ARC  BODY  CURVE.  WHEN  KSI  *  1.  THE  BODY' 
420  PRINT  TAB(20> ► 'CURVE  IS  THE  ARC  OF  LEAST  RADIUS  OF  CURVATURE' 

430  PRINT  TAB(20) » 'WHICH  IS  TANGENT  TO  THE  BODY  AT  X=0' 

600  FOR  1=1 T05 
610  PRINT 
620  NEXT  I 

625  PRINT  'DEFAULT  VALUES  OF  INPUT  VARIABLES' 

626  PRINT 

630  FOR  1=1 T05 

640  PRINT  *  *»A*(I). 

650  NEXT  I 
660  PRINT 

670  PRINT  R1.R2.L.P.G9 
680  PRINT 

690  FOR  I  =  6  TO  10 
700  PRINT  '  ' ? A$ ( I > . 

705  NEXT  I 
710  PRINT 
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720  PRINT  K5* A0*M*N*N1 
730  PRINT 

740  PRINT  *  *  '  *  *A$< 11 ) 

750  PRINT  K9 
800  PRINT 

810  PRINT  'NEW  VARIABLE? 

820  INPUT  B**Z1 

830  FOR  1=1  TO  13 

850  IF  A*<I)  =  B$  60  TO  890 

860  NEXT  I 

870  PRINT  'VARIABLE  NOT  RECOGNIZED*  TRY  AGAIN' 

880  GO  TO  810 

890  IF  B»=A*(12)  GO  TO  810 

900  IF  B*  =  A$< 13)  GO  TO  1160 

910  ON  I  GO  TO  920* 940*960* 980* 1000* 1020* 1040* 1060* 1080? 1100* 1120 

920  R1=Z1 

930  GO  TO  1150 

940  R2=Z1 

950  GO  TO  1150 

960  L=Z1 

970  GO  TO  1150 

980  P=Z1 

990  GO  TO  1150 

1000  G9=Z1 

1010  GO  TO  1150 

1020  K5  =  Z1 

1030  GO  TO  1150 

1040  A0=Z1 

1050  GO  TO  1150 

1060  M=Z1 

1070  GO  TO  1150 

1080  N-=Z1 

1090  GO  TO  1150 

1100  N1=Z1 

1110  GO  TO  1150 

1120  K9=Z1 

1150  GO  TO  810 

1160  J1*SQR((P"( <G9-1)/G9)-1)*2/<G9-1) ) 

1161  IF  J1>1  GO  TO  1170 

1162  T =0 

1163  R3=R2 

1164  GO  TO  1190 

1170  T=SQR( <G9+1)/<G9-1) )*ATN<SOR( <G9-1 )/ ( G9+1 > * ( Jl* Jl-1 ) ) ) -ATN ( SGR ( J1*J1-1 ) ) 
1180  R3=R2*P"<  <G9+l)/2/G9)*(2/(G9+l >  >" ( ( G9+1 ) 72/ ( G9- 1 ) ) / Jl 
1190  J2=A0*J1/SGR<1+<G9-1)/2#J1*J1> 

1200  R9=<Rl-R2+L*L/<Rl-R2>>/2 
1210  R9=R9*R9 

1220  L.9=SQR (  ( R1-R2 ) #<R1-R2 ) +L*L  ) /2 
1230  S9=SQR(R9)-SQR<R9-L9*L9) 

1235  IF  K9=0  GO  TO  1285 
1240  S=S9*N9 
1250  R8=<S*S+L9*L9>/2/S 
1260  G1=ATN(L/(R1-R2) ) 

1270  Hl=L/2-<R8-S)#C0S<Gl ) 

1280  Kl=(Rl+R2)/2-<R8-S)#SIN<Gl) 

1282  IF  Jl  >  1  GO  TO  1285 

1283  R6=123456 
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1284  GO  TO  1290 

1285  R6=(Rl-R2)/< l-COS(T) ) 

1290  H2=L+R6#SIN<T> 

1310  K2=R3-R6 
1320  B2=1-M*M 
1330  U9=1117*M 
1340  J5=J2/U9 
1350  X1=0 

1355  C1=0 

1356  01=0 

1360  GOSUB  3000 
1370  01=R 

1500  PRINT  'INPUT  VARIABLES  FOR  THIS  RUN* 

1510  PRINT 

1520  FOR  1=  1  TO  5 

1530  PRINT  '  afA«<I>» 

1540  NEXT  I 

1550  PRINT  R1 fR2»L»P»G9 
1560  PRINT 

1570  FOR  I  =  6  TO  10 
1580  PRINT  '  *  $ A$< I ) > 

1590  NEXT  I 

1591  PRINT* 

1600  PRINT  K5»A0»M»N»N1 
1610  PRINT 

1620  PRINT  *  " i At < 1 1 > 

1630  PRINT  K9 
1640  PRINT 

1650  PRINT  'CALCULATED  PARAMETERS  FOR  THIS  RUN' 

1660  PRINT 

1670  PRINT  *  R-P'»*  THETA-J*  t '  BODY  RADIUS' r*  PLUME  RADIUS* f'  PLUME  TERr 

1680  PRINT  R3fTfR8fR6fH2 
1690  PRINT 
1700  PRINT 

1710  PRINT  *  X'f'  R * »  *  U/U-O'f*  V/U-O'f*  C-P'f'  C-D  INT 

1720  PRINT 
1730  R8=R8*R8 
1740  R6=R6*R6 
1750  M1=(R2-R1)/L 
1790  D1=L/N 
1795  I9=N/10 
1800  FOR  I  =  0  TO  N 
1805  X=D1#I 
1810  11=0 
1820  X8=L/N1 
1830  12=0 
1840  X1=X 
1850  GOSUB  3000 
1860W3=R*R*B2 
1870  X1=0 

1880  GOSUB  3000 

1881  GOSUB  3500 

1882  I1=I1+F 
1884  I2=I2+G 
1894  X8=L/N1 

1896  IF  H2<10*L  GO  TO  1902 
1898  I5=10*N1 
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1900  GO  TO  1920 

1902  I5=H2/X8 

1903  IF  JlOl  THEN  I5=N1 
1920  FOR  13=1  TO  15-1 
1925  X1=I3*X8 

1930  GOSUB  3000 
1935  GOSUB  3500 
1940  14=13/2 

1945  IF  <14-1 NT (14))/. 5> .5  GO  TO  1965 

1950  I1=I1+2*F 

1955  I2=I2+2*G 

1960  GO  TO  2010 

1965  I1=I1+4*F 

1970  I2=I2+4*G 

2010  NEXT  13 

2011  X1=(I3-1)#X8 

2012  GOSUB  3000 

2013  GOSUB  3500 

2014  I 1=I1+F 

2015  I2=I2+G 

2016  11=1 l#X8/3 
2018I2=I2*X8/3 
2020  X1=X 

2030  GOSUB  3000 

2040  Il=ll/2+K5*R2/4*<J5  1 )/SQR(X*X+B2*R*R> 

2050  I2=I2/2*B2*R-K5/4*R2/R* ( J5-1 > * ( X/SGR ( X*X+B2*R*R ) +1 ) 
2060  C=-( 2# 11+11*11+1 2# 12) 

2090  S7=R-Q1 

2100  C1=C1+2*3.14159*R*S8*D1*C 

2190  C3=Cl/3 . 14159/R1/R1 

2195  IF  I/I9-INT<  1/19)00  GO  TO  2210 

2200  PRINT  XrRf Il»I2fC»C3 

2210  Q1=R 

2220  NEXT  I 

2230  FOR  1=1  TO  5 

2240  PRINT 

2250  NEXT  I 

2260  PRINT  "DO  YOU  WANT  A  NEW  RUN?  TYPE  YES  OR  NO. 

2270  INPUT  B* 

2280  IF  B*=A*( 12)  GO  TO  810 
2290  STOP 

3000  IF  XI >L  GO  TO  3020 
3005  IF  K9=0  GO  TO  3040 
3010  D2=H1 -X 1 
3012  D3=D2*D2 
3014  D4=SQR<R8-D3) 

3016  R=K1+D4 

3018  S8=D2/D4 

3019  RETURN 

3020  D2=H2-X1 
3022  D3=D2*D2 
3024  D4=S«R<R6-D3) 

3026  R=K2+D4 

3028  S8=D2/D4 
3030  RETURN 
3040  R=M1*X1+R1 
3045  S8=M1 
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3050  RETURN 
3500  W=X~X1 
3510  W1=W*W 
3520  W5=SQR<W1+W3) 

3530  W4-R*SP/W5/< W1+W3) 
3540  F=W4*W 
3550  G=W4 
3560  RETURN 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 


1.  DTNSRDC  REPORTS.  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH¬ 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS,  A  SEMIFORMAL  SERIES,  CONTAIN  INFORMATION  OF  A  PRELIM¬ 
INARY,  TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES.  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN¬ 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE-BY-CASE 
BASIS. 


